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Abstract 

Flat directions charged under an R-symmetry are a generic feature of O'Raifeartaigh 
models. Non-topological solitons associated with this symmetry, R-balls, are likely to 
form through the fragmentation of a condensate, itself created by soft terms induced 
during inflation. In gravity mediated SUSY breaking R-balls decay to gravitinos, 
reheating the universe. For gauge mediation R-balls can provide a good dark matter 
candidate. Alternatively they can decay, either reheating or cooling the universe. 
Conserved R-symmetry permits decay to gravitinos or gauginos, whereas sponta- 
neously broken R-symmetry results in decay to visible sector gauge bosons. 



1 Introduction and conclusions 



Of upmost importance in any supersymmetric theory is the appearance of classical flat 
directions, or pseudo-moduli: scalar combinations along which the classical potential is 
invariant. Indeed, ref. [T] demonstrates that the quantum mechanical stabilisation (or not) 
of flat directions is instrumental in determining whether a SUSY breaking vacuum is viable. 
Furthermore any renormalisable O'Raifeartaigh model [2] necessarily has such a direction, 
namely the scalar component of the goldstino superfield [3111]. 

Of course, flat directions are well known to exist in the MSSM [5] and are frequently 
utilised when explaining the source of the universe's baryon asymmetry [5J. High scale 
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SUSY breaking associated with inflation can result in a condensate with large expectation 
value forming along a flat direction charged under baryon number. If the model does not 
preserve it at high energies, this condensate carries significant fractional baryon number 
that survives to the present day. Originally it was thought that the condensate would evap- 
orate into light fermions (imparting its baryon number to visible sector degrees of freedom) 
but later studies [7HH] instead suggest that the condensate fragments into localised lumps 
known as Q-balls [TB] . 

Given that O'Raifeartaigh models have flat directions of their own, one can ask whether 
something similar happens in the hidden sector. Recent work [T6H20] has certainly sug- 
gested that other, non-MSSM flat directions can have important consequences. Nelson 
and Seiberg's powerful argument [21] implies that a generic, calculable model can only 
break M = 1 SUSY if it has a (possibly spontaneously broken) U(l) R-symmetry. Even if 
the SUSY breaking is only metastable one expects an approximate R-symmetry, with any 
R-violating operators being suppressed to ensure a sufficiently long lived vacuum. Rather 
than baryon number, it is this R-charge that is carried by flat directions in O'Raifeartaigh 
models, thus leading to R-balls: non-topological solitons stabilised by a global U(l) R- 
symmetrjQ. 

In this paper we investigate the general, cosmological evolution of flat directions in 
O'Raifeartaigh models. Since many such models of SUSY breaking are only low energy 
effective descriptions, it is important to consider effects from the underlying microscopic 
theory. The cutoff scale of the low energy theory will prove to be an important parameter. 
Initially, the flat directions evolve similarly to their MSSM counterparts and a spatially con- 
stant condensate is formed. As the universe cools down the condensate fragments (subject 
to constraints imposed by thermal effects) and forms R-balls. At this stage it doesn't make 
a difference whether R-symmetry is spontaneously broken or not. The effective potential 
for the flat direction preserves an R-symmetry over the majority of its domain and only 
small field values, i.e. the edges of the R-ball, see the effects of any breaking. Eventually 
R-balls decay, either by surface evaporation to light fermions if R-symmetry is preserved, 
or to light bosons throughout the interior if R-symmetry is spontaneously broken. 

Two distinct phenomena can then be attributed to R-balls. If sufficiently long lived 
they act as cold dark matter and are easily capable of providing the dominant contribution 
to the observed density. This scenario only occurs in gauge mediated SUSY breaking where 
the SUSY breaking scale (or equivalently the gravitino mass) is small enough for R-balls 
to outlast the current age of the universe. Otherwise R-balls decay. A significant density 
during decay leads to a second bout of reheating or a cooling of the universe, depending on 
the temperature of the decay products. Allowing the universe to pass through an R-ball 
dominated epoch thus decouples the dynamics of the inflaton from the generation of visible 
sector matter. Gravity mediated SUSY breaking sees R-balls evaporate to gravitinos, re- 
heating the universe. Gauge mediation allows for both reheating and cooling, by gravitinos 

lr This may not be the only bearing global symmetries have on the cosmology of the SUSY breaking 
sector. Ref. [22] entertains the possibility of soliton induced, SUSY restoring phase transitions in models 
of metastable SUSY breaking. 
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and gauginos if R-symmetric, and by visible sector gauge bosons when R-breaking. The 
parameter space of some example O'Raifeartaigh models is dissected in figures IBTTj 15 . 2 \ I5~3"l 
In each case there remain several regions where R-balls are forbidden, either by thermal 
effects or due to contradiction with experimental evidence, or form but have little obvious 
effect. 

The remainder of the paper is organised as follows. In section [2] we will review the 
observation of refs. [2111] regarding the existence of flat directions in O'Raifeartaigh mod- 
els, and will discuss how they are lifted by loop corrections. Section [3] deals with the 
cosmological evolution of the flat direction, from inflation to fragmentation, including po- 
tential thermal effects. R-balls themselves are investigated in section |4j with a particular 
emphasis on their decay rates. Finally, in section [51 we will consider the phenomenological 
consequences. Throughout this work an attempt is made to keep the analysis as general 
as possible. It is inevitable that certain aspects (decay, for example) would benefit from 
a more specific, numerical treatment. However, the aim here is to broadly consider the 
entire spectrum of consequences associated with flat directions in O'Raifeartaigh models. 



2 Flat directions in O'Raifeartaigh models 

To begin we will recap the results of refs. j3jH] on flat directions in O'Raifeartaigh models, 
followed by a discussion on their lifting. Consider a model of SUSY breaking which is 
generic, calculable and has a low energy description in the form of an O'Raifeartaigh 
model. Ref. [21] tells us that this model must have an R-symmetry. Denoting the chiral 
superfields as {$ a } the scalar potential is given by the F-terms 

Z7($) = FlF a = WlW a (2.1) 

where W a is the derivative dW/d<f> a . For a locally stable SUSY breaking vacuum to exist 
we must be able to find a solution $ a = ($) a to the equations 

dU 

__ = W ab W ] b =0 iy a ^0 for at least one value of a. (2.2) 

C» a 

The first equation contains the fermion mass matrix M.1/2 = W ab so these two statements 
confirm that there is at least one massless fermion, the goldstino, corresponding to the 
fermionic component of the chiral superfield whose F-term gets a VEV. Meanwhile the 
scalar mass matrix is 

Ml = ( M\ /2 M lft m ^ j where F ^ = M (2 3) 

and must be positive semi-definite for $ a = ($) a so that the vacuum is free of tachyons. 

Now take the direction X a = WJ|$ 0= ($) (the scalar component of the goldstino super- 
field) and consider the norm 



(2.4) 
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This quantity necessarily vanishes. Otherwise one could always rotate the phase of X to 
make the right hand side negative, contradicting the assumption that the vacuum is locally 
stable. Ergo it must be that TX = 0, therefore (X, X*) is a massless scalar. In fact the 
claim is much stronger than this for a renormalisable superpotential. TX = can be 
expanded as 

(TX) a = \w abc Wlwl\ = 0. (2.5) 

Moving along the X direction the scalar potential changes as 

[/«$) + AX) = [/(($)) + 5W\5W a (2.6) 
for some complex parameter A, where 



5W, 



a 



&W ab Wl + l -A 2 W abc Wlwl 



(2.7) 



This expansion is exact up to fourth derivatives in W (i.e. for a renormalisable superpo- 
tential) and, upon consulting eqs. (12. 2 p and (12. 5p . one immediately finds 5W a = 0. In 
other words the scalar potential is invariant under complex translations along the scalar 
component of the goldstino superfield: it is a flat direction. For non-renormalisable models 
the scalar partner of the goldstino remains massless, but the degeneracy may be lifted by 
higher order terms in the superpotential. 

A consequence is that any O'Raifeartaigh model can be recast in terms of a goldstino 
superfield X, whose F-term VEV is responsible for SUSY breaking, and some other super- 
fields {(fi} [I]. In this basis we will allow the most general, renormalisable form for the 
superpotential 

W = fX + (Hij + \ ij X)Lf i Lp j + K ijk LpiLpjLp k (2.8) 

for coupling constants /, /i, A and k, where / is assumed real and positive without loss of 
generality. Note that these symbols will frequently be used without their indices to denote 
the generic size of the couplings. The superfields <fi are defined so as to have vanishing 
VEVs whereas the scalar component of X is our classical flat direction. Transforming into 
this basis may not respect the global symmetry group of the model but, in a model with 
an R-symmetric vacuum, X clearly has R-charge +2. Any renormalisable O'Raifeartaigh 
model with an R-symmetric vacuum therefore possesses a flat direction with non-zero R- 
charge. 



2.1 Lifting the flat direction 

For the vacuum to be well defined X cannot remain flat and must be stabilised. At low 
energy this is accomplished by quantum effects, typically at one loop via the Coleman- 
Weinberg potential [23] evaluated with respect to a UV cutoff scale M. Around this scale 
we may start to see non-renormalisable effects from an underlying microscopic theory, 
which could also lift the flat direction, but are highly suppressed at low energy so have 



4 



little impact. Exactly where the Coleman- Weinberg potential is minimised depends on the 
details of the model but one can deduce its approximate form in two limits. Close to its 
minimum at X = (X) the effective potential goes like 



U eS (X) = U + -m 2 \X - (X)\ 2 + Q(\X - (X)\ 3 ). 



(2.9) 



If R-symmetry is preserved this minimum must be at (X) = 0. The mass term m 2 depends 
on the couplings /, /i and A and one can deduce its value either by explicitly calculating 
the Coleman- Weinberg potential for a given model, or by using the more elegant methods 
developed in ref. |24) 2 I . Schematically it should go like 

A 4 f 2 

m ' ~ vd? < 2 - 10 > 

At the other end of the scale when X is large, specifically \fi + XX \ 2 ^> Xf, the SUSY 
breaking is small. In fact an appeal to naturalness suggests / and fi 2 should not be too 
dissimilar so small SUSY breaking often implies \XX\ ^> \fj,\. In this limit we can find an 
alternative formulation for the effective potential by integrating out the massive tp's. The 
process contributes an extra term to the Kahler potential 

K eS (X,X^) 



\X\ Z - 



Tr 



and consequently leads to a one loop effective potential 



U eS (X) 



x 



d x d x \K eS 



16tt 2 



5> 



1 + ln 



\kX\ 
M 



0(X 4 



(2.11) 



(2.12) 



The mass matrix Ai = XX has been taken from the superpotential (12. 8p (neglecting /x 
in the limit of small SUSY breaking) where {Aj} are the eigenvalues of A. If fi 2 3> / for 
some reason, the Kahler potential is modified by the possible appearance of additional 
light states at certain points in the pseudo-moduli space, which cannot be integrated out. 
These lead to singular behaviour and our approximation must be reassessed. Regardless, 
all conclusions reached apply wherever <C \XX\ <ti M. Note also that eq. (I2.12p is 
only valid up to second order in / and A. If the A's are large, or there are a vast number 
of non-goldstino superfields, higher loop corrections become important and invalidate the 
Coleman- Weinberg formula. 

In the following we will be interested only in the potential relative to the SUSY breaking 
vacuum so will omit the constant term of eq. (12.91) . which goes like f 2 up to loop effects. 
Normalising both potentials thus gives the final form 

e(X) w -m 2 |Y- (X)| 2 forX 2 </ 



U, 



U eS (X) 



167T 



C + hi 



\XjX\ 
M 



for / < X 2 < M 2 



(2.13) 



2 It was also shown in ref. 
not equal to or +2. 



that spontaneous R-symmetry breaking requires a superficial with R-chargc 
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C/eff 




Figure 2.1: A typical low energy effective potential for the scalar component of the goldstino 
superfield X. At tree level X is a flat direction but is lifted by loop effects. For X 2 <C f 
these are polynomial in nature and stabilise X at some value (X) (equal to zero in an 
R-symmetric vacuum). For X 2 3> / the quantum corrections become logarithmic. 



for some order one parameter (. For a model that is perturbatively well behaved and whose 
couplings satisfy / ~ /i 2 this effective potential is illustrated in figure 12. 1L 

3 Cosmological evolution 

Flat directions in the MSSM have been extensively studied (see e.g. ref. [5] and the refer- 
ences therein) and many of the techniques used are directly applicable here. In particular, 
a condensate forms along flat directions in a wide variety of models due to the Affleck-Dine 
mechanism [B|. The idea is that quantum fluctuations are spread out during an inflationary 
period of the universe's evolution, with only the long wavelength modes surviving to form 
a spatially constant condensate. Initially, the expectation value of the condensate field is 
set by high scale, SUSY breaking couplings to the field driving inflation: the inflaton. As 
the universe cools down, the expectation value of the condensate decreases until inflaton 
effects become subdominant and the condensate field begins to move in the low energy 
effective potential. At this point the condensate becomes unstable to spatial perturbations 
and begins to fragment [THTij. the final state being non-topological solitons stabilised by 
some conserved charge: Q-balls [T5] . 

In O'Raifeartaigh models the overall picture is similar. We will keep with standard 
nomenclature and refer to the flat direction X as the condensate field. A summary of the 
key points is then as follows. 

• Inflation: A tachyonic soft mass, originating from the high scale SUSY breaking 
driving inflation, can drive the condensate field away from its low energy VEV. The 
condensate field is stabilised by higher order soft terms originating from a microscopic 
theory. For an accidental R-symmetry the condensate is stabilised at a scale A, 
parametrically between M and y/f, whereas for an exact R-symmetry it is stabilised 
around the Planck scale. 
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• Rotation: The Hubble parameter decreases until the inflaton induced soft terms are 
comparable in scale to the low energy effective potential. The low energy vacuum is 
restored and the condensate field begins to rotate in its potential well with magnitude 
A and frequency Xf/AnA. 

• Fragmentation: Shortly after beginning rotation the condensate is rendered unsta- 
ble to spatial perturbations and fragments into localised lumps of size 

The fragments coalesce into extended, classical objects with large charge: R-balls. 

• Decay: In R-symmetric models R-balls evaporate to gravitinos at tree level and 
other light fermions at one loop. In models with spontaneously broken R-symmetry 
R-balls can decay more quickly into other light fermions at tree level and also into 
light bosons. 

We shall now discuss each stage of the evolution in more detail. 
3.1 Inflation 

At such early epochs one cannot consider the low energy O'Raifeartaigh model in isolation. 
Inflation induces extra soft terms in the potential related to the Hubble parameter |26j . 
Non-renormalisable operators arising from the microscopic theory also have a profound 
effect on the early universe dynamics, especially if they do not respect the R-symmetry of 
the low energy theory. 

Cosmologically the scalar component of X obeys the equation of motion 

X + 3HX - \\/ 2 X + ^ = (3.1) 

or oa t 

for a high energy effective potential V e s(X). H is the Hubble parameter, a is the scale 
factor of the universe and, for a homogeneous condensate, the gradient term obviously 
disappears. In the absence of thermal effects (the discussion of which is postponed until 
section l3~4"l) the general form of V e g(X) is known [26] to be 

V cS (X) = -cH*\X\> + j^r 3 (Ar,X n + h.c.) + jJ^lvlW^ (3.2) 

for a microscopic superpotential coupling 77, order one constants^) c and A, and where n > 4. 
The first term is a soft mass induced through inflaton couplings and is always present. The 
remaining terms are generated by non-renormalisable superpotential operators originating 
from the microscopic theory, which has been allowed to break R-symmetry (i.e. the R- 
symmetry of the low energy theory is accidental) and lift the flat direction. Even if R- 
symmetry is exact up to the Planck scale it is expected to be broken by gravitational 

3 For D-term inflation one expects ^4 = 0. This has little effect for an accidental R-symmetry but, 
when R-symmetry is exact, D-term inflation precludes the possibility of R- violating operators and the 
condensate always has vanishing charge. 
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effects. If no non-renormalisable superpotential operators lift the flat direction, one would 
instead have 2 

VMX) = -cH 2 \X\ 2 + -^il"- 1 + -f-^5|X| 2n - 4 (3.3) 

ettV ) II M n-3 M 2n-2 II V J 

but the conclusions of the subsequent discussion are unaffected, other than replacing the 
scale A (defined shortly) with the Planck scale Mp. 

The soft mass in eq. (13. 2p is extremely important in the early universe when H is large. 
For the minimal Kahler potential K = X<X it is generated by supergravity corrections of 
the form 

V eS (X) = e K ^V ini ( X ) (3-4) 

where \ is the inflaton superfield. During inflation the inflaton vacuum energy dominates 
the universe so VJ n f (x) ~ H 2 Mp, leading to a soft mass term with negative c. Therefore the 
effective potential is minimised at the origin and the dynamics are uninteresting. However, 
for a non-minimal Kahler potential it is quite possible that c is positive. Consider, for 
example, the term K D (x^X){X^X)/Mp which is allowed by all possible symmetries of 
the model. In fact terms of this form are inevitable for superpotentials like (\2.8\i : they 
arise as counterterms for Yukawa couplings [261428] . Their contribution to the soft mass is 



J M p M p 



(3.5) 



so a positive Kahler potential coefficient results in a negative coefficient for the soft mass 
term, potentially winning out over the previous contribution. Since this soft mass is gen- 
erated by the high scale SUSY breaking associated with inflation rather than the SUSY 
breaking associated with the MSSM, the result is independent of the scale appearing in 
the O'Raifeartaigh model itself. 

Assuming the soft mass term is tachyonic, the effective potential is initially unstable 
around the origin and is only stabilised at 

X ~ (HM n - 3 ) 1/(n - 2) V cS (X) ~ (#n-l M -3)2/(n-2) ^ 

with a choice of n — 2 distinct minima corresponding to different choices of phase. The 
condensate field quickly settles into one of these minima and remains there throughout 
inflation due to the large, Hubble induced damping term in eq. (13.11) [2"E] . Immediately 
after inflation the Hubble parameter evolves as 1/1 The minimum thus moves closer to the 
origin over time until V e s(X) ~ A 2 / 2 /167r 2 . At this stage the low energy effective potential 
(|2~T3|) takes ovefl the corresponding Hubble parameter and condensate expectation value 



4 We shall assume there are no independent minima at large X, i.e. ones that do not require a negative 
mass term centred around the origin to be stable. If such minima did exist, the condensate expectation 
value could remain near the cutoff scale and the theory would never flow into its low energy O'Raifeartaigh 
model description. This effect is not of interest here, but could provide a novel mechanism for models with 
uplifted vacua to find themselves in a higher energy vacuum. 
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being 



H 



1 f^P 



M n ~ 3 \ 4vr 



I/< "" 1 ' - ^A,^)^". (3.7, 



Here, we have defined the parameter A which will be important in all that follows. This 
scale is parametrically between the cutoff scale of the O'Raifeartaigh model M and the 
loop suppressed SUSY breaking scale \J Xf /Air. It actually turns out that the case A 2 < f 
is uninteresting so we will henceforth assume / 2 < A < M. It will also be convenient to 
recast the above value of the Hubble parameter in terms of A 

H ~T^- ( 3 - 8 ) 
47rA v ' 

3.2 Rotation 

Below H ~ Xf /AirA the condensate performs rotations about the minimum of the low 
energy effective potential (I2.13p . Whether R-symmetry is spontaneously broken or not, the 
logarithmic regime of eq. (12.13]) is independent of the condensate's phase. An effective R- 
charge is thus conserved for the most part, resulting in approximately circular or elliptical 
motion. The charge stored in the condensate is determined by the interplay between 
the various R-violating operators in the high energy effective potential (13. 2p when the 
rotation begins. At the transition point between high and low energy regimes these are 
of comparable size to the R-preserving terms so can impart a sizeable 'torque' on the 
condensate and bestow it with a large fractional charge^). The only subsequent source of 
R-symmetry violation is the small X regime, hence we can think of the trajectory as being 
smooth and elliptical, but possibly getting a kick if it gets too close to the bottom of the 
potential well. 

At first X 2 ~ A 2 3> / so the large X limit of eq. (I2.13P applies. At any given time 
one could ignore the damping term in eq. ( 13. ip and find a circular solution X = X c e lut , for 
constant amplitude X c and frequency 

2 _ J_dUeS X 2 f 2 

Elliptical variants will be mentioned in section I4TT1 Initially X c ~ A and H ~ Xf/AnA so 
the damping coefficient is similar to the frequency and the motion is critically damped. As 
the Hubble parameter continues to decrease, the damping follows suit and the motion be- 
comes underdamped. We thus expect the effect of the damping to be small, the condensate 
expectation value and frequency of rotation remaining around A and Xf/AwA respectively. 



5 Relating R-charge to baryon number could thus lead to some interesting asymmetric dark matter 
scenarios such as those in refs. [T5ll20j . 
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3.3 Fragmentation 



Left to its own devices the damping would, eventually, force the condensate into the small 
X regime of eq. (I2.13P and it would continue to perform underdamped oscillations about 
X = (X) with frequency m and magnitude yff. However, over the domain / <C X 2 < A 2 
the effective potential is logarithmic so increases slower than quadratically. A condensate 
oscillating in this kind of potential behaves as matter with a negative pressure, i.e. it is 
unstable with respect to spatial perturbations [7HH]. 

For circular rotations it is possible to see this explicitly and estimate the typical size of 
the fragments [EJ[T2]. Working in the underdamped regime discussed above one can write 
down the approximate solution 

X w A e ^//4-A)^ (3 10 ) 

Now consider fluctuations in the magnitude and phase of X, of the forms <5£ = 8£oe at+lkx 
and 56 = 56oe at+lkx respectively. Unstable modes have a > 0. Substituting into the 
equations of motion a non-trivial solution for <5£ and 68q exists only if 

a + 2 77 + 7777777 « + 77 " 777777 ) k = C 3 - 11 ) 



a 2 16tt 2 AV V^ 2 8n 2 A 2 
which is satisfied for real, positive a and k only if 

'k 2 \ 2 f 2 \ n 



< < - < (3.12) 

2 8tt 2 A 2 ; a 2ttx/2A 1 ; 

The most amplified mode (largest a) has a scale set by the upper limit of this inequality, 
so the fragments have typical size 

a 27ry/2A 

r~-~^-. (3.13) 

Moving away from circular trajectories requires a numerical approach, but for a logarithmic 
effective potential it only changes the result by factors of order unity [T2]. Regardless, 
fragmentation occurs after the horizon size if -1 surpasses the typical size of the fragments. 
This has already happened when the condensate begins to rotate at H ~ Xf /An A so the 
process can take place immediately. 



3.4 Thermal effects 

So far we have ignored the fact that all of these processes take place in a thermal bath. 
Below the reheat temperature Th the inflaton decays into light degrees of freedom. This 
reheats the universe, at which point it can be considered a cooling plasma with temper- 
ature T ~ y/HMp. Even beforehand the universe contains a dilute plasma with a lower 
temperature T ~ (HMpT%) lj/4: |37ffi The permitted values for the reheat temperature are 

6 Note that this temperature should not be associated with the temperature of the R-balls as they are 
not generally in thermal equilibrium. 
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constrained by gravitino cosmology. If reheating occurs too early gravitino decay prod- 
ucts disrupt nucleosynthesis or, if the gravitino is stable, its relic density may be too 
high [38144"2] . For an unstable gravitino (e.g. gravity mediation) one requires Th < 10 6 m3/2 
whereas for a stable gravitino (e.g. gauge mediation) the bound is Th < 10 14 (m3/2/GeV) -2 
GeV. Additionally, there is a lower bound of 10 MeV placed on the reheat temperature to 
ensure that reheating happens before nucleosynthesis. 

Any particle coupling to the condensate field at tree level acquires a mass of order | XX | . 
If this mass is greater than the plasma temperature, i.e. |AX| > T, the state decouples 
from the condensate. Otherwise the condensate is eroded. The condensate expectation 
value goes like A during fragmentation so one must have 

A > ^ (3.14) 

for the process to proceed unhindered, where Tf is the temperature of the universe at this 
time. If the inequality is satisfied there remain thermally sensitive loop corrections to the 
the low energy effective potential (I2.13P to take into account [T2]. Couplings to heavy 
states yield the extra term0 

\2t^4 /\x\ 2 \ 
U cS (X,T) = — -In . (3.15) 



16vr 2 V T 2 , 

Unless f > T 2 this contribution dominates and the above expression should be used in 

place of eq. ( I2.13p . However, since the form of the effective potential (which arises from 
loop corrections anyway) is not changed by the thermal contribution, its effects can be 
absorbed into the parameter / 

[ f for T 2 < f 

f — > h = < (3.16) 

J It 2 for t 2 > f. 1 ; 

Other thermal effects exist but are expected to be small for the large R-balls we will be 
discussing [8JII3] so will not be considered. 

According to section 13.31 the condensate fragments at H ~ A/t/47tA so, using the 
expressions given at the start of this section to relate the Hubble parameter to the plasma 
temperature, this equation breaks up into four possible domains depending on the relative 
sizes of the scales 



D 3 



T f <T h ,T 2 <f Tj/M P ~ A//4ttA 

T f >T h ,T]<f T^MpTl ~ A//4vrA 

Tf > T h , Tf > f =► Tf/MpTj 2 ~ XT] /An A 
T f <T h ,Tf > f Tf/Mp ~ XT]/AttA. (3.17) 



7 In addition one should include a temperature dependence in A. We shall assume this is small enough 
to omit here so as to avoid further complicating the discussion. 
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There are then three different solutions for the fragmentation temperature; the fourth 
solution is just the boundary between two other domains. In each case the important 
parameters can be summarised as follows 



Domain 



D 3 : 



A > max 



An 



MMp > A > XMpT h 



4vrT, 



A < min 



Airf 

XMpT* XM P 

Anf ' An 



T 



I 



P. 



T 



V 47rA 

AnA J 

XMpTlV' 2 
AnA J 



f 
f 



\M P Tl 
AnA 



A, 



( fM P \ 
\ AnX ) 



TJT 



4ttA 3 
\ AnX 



(3.18) 



where A min is the value demanded by the inequality f)3.14p . 



4 R-balls 

Much as MSSM condensates end up as B-balls (or L-balls), our condensate collapses into 
extended, classical configurations with large charge, R-balls, which are non-topological 
solitons formed under the influence of the U(l) R-symmetry [TJ)J[22J|31]. Recall that the 
low energy effective potential (12. 13[) preserves R-symmetry in the large X regime whether 
or not it is spontaneously broken in the vacuum, hence R-balls form in both cases but have 
different decay properties depending on the eventual status of the symmetry. 

R-ball solutions take the form X(x,t) = X(x)e luJt , for a real parameter u and a real 
function X(x) that minimises 

d 3 x (\\WX\ 2 + U U1 {X) j (4.1) 



where a new potential 



2' 



U U (X) = U eS (X) - Ilo 2 \X\ 2 (4.2) 



2 

has been defined. Assuming spherical symmetry the associated equations of motion are 

d ! X + 2dX_d«, = 
dr 2 r dr dX 

for a radial coordinate r. One can always find a solution if U UJ (X) meets two criteria: it 
retains a stable minimum at X = (X) and there exists a non-zero X such that U^X) < 
U U ({X)). As long as oj < m the potential defined by eq. (I2.13P satisfies both; uj 2 X 2 always 
beats the logarithm for large X. 

For the case at hand we can be more precise owing to the effective potential being so 
flat [151I1B]. Substituting eq. (12.131) into the equation of motion (14.31) yields two limits: 

9 f m 2 X for X 2 < / 

X" + -X' + u 2 X = I (4.4) 
r ] CA 2 /|/16vr 2 X for X 2 > /. V 1 
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Figure 4.1: A typical R-ball solution. At small and large r it has the forms sin ur/ur 
and e~ mr respectively, with a width of about tt/uj. The parameter uj scales like _R _1//4 , the 
height like i? 1//4 and the energy i? 3//4 . 



In the first, the equation is exactly soluble. In the second, the right hand side can be 
ignored as X is large, in which case the equation is again exactly soluble. Choosing the 
solution that is finite at r = 0, decays to zero as r — > oo and is continuous between the 
two regimes we find 

X = — sin cur for r < , X = — sin ur R e^-OVm*-^ for r > Tr _ 
ujr uir 

The R-ball's width is approximated by r R m tt/uj, such that the cross over occurs near 
the first zero of sinwr where X is arbitrarily small. Note that a cross over to large X is 
inevitable because the small X solution diverges as r — > 0. When u m the R-ball is 
very large and the solution can be further simplified to find 

X = sinwr for r < r R , X = sinwr^ e m{ - rR ~ r> for r > tr. (4.6) 

cor tur R 

Evaluating the energy of this field configuration and choosing ur R « ix gives 

E = + ouR + 0(n - ur R ). (4.7) 

Minimising with respect to u and explicitly evaluating the charge finally yields expressions 
for the R-ball parameters 



u R ~R-V 4 ,/Xfr r R ~—= XR^R^y/Xfr E R ~R 3 /*y/Xfr (4.8) 

V^JT 

up to order one coefficients. Putting it all together, a typical R-ball solution is illustrated 
in figure 14.11 
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4.1 Formation 



In section 13731 it was shown that, for an initially circular trajectory, the condensate forms 
fragments of size r ~ 27ta/2A/A/t shortly after it starts to rotate with magnitude X c ~ A 
and frequency v c ~ A/t/47tA. We can use these details to estimate the typical charge [12]. 
The charge density during fragmentation is given by 2v c Xl so the total charge per fragment 
is 

100A 4 tl . 

R-^jf- (4-9) 

If each fragment collapsed into a single R-ball this would be their charge. Actually, numer- 
ical simulations [12] suggest that R-balls form slightly after fragmentation in a logarithmic 
potential, leading to a reduced charge of 



Each fragment therefore contains about 10 3 R-balls. 

Using the energy (14.81) the energy density stored in R-balls just after formation, at 
H ~ A/r/47rA, can now be determined to be 

p R ( Tf )^^Il =^ Q R (T f ) = pR{tf) ~ — (4 11) 

PR{ f) 10 R{ f) 3H*M 2 P M 2 P 1 } 

unless A ~ Mp whereupon Qp ~ 1 and R-balls dominate the universe until they decay. 
Of course, this result scales with the expansion of the universe. If R-balls form above a 
temperature of y/f there is also lowering in density due to the decreasing value of fr (I3.16p . 
Treating R-balls as non-relativistic matter the density can be related to the scale factor 
via pR ~ a -3 , implying that 



Jfr for matter domination 

ftii(T) oc { (4.12) 
s/Jt/T for radiation domination. 

It should be noted that there are no other processes to form R-balls once the condensate has 
fragmented so they cannot maintain thermal equilibirum. The idea of solitosynthesis [4T1 - 
4"5] does not apply here as the goldstino carries R-charge and is lighter than the condensate 
field. 

To generalise to more eccentric condensate trajectories one can define the parameter 
e — vjv c — knvhj Xfx'- the ratio of the angular velocity of the condensate to its maximal, 
circular value. Equivalently it can be thought of as the fraction of the maximum possible 
charge that is stored in the condensate. Clearly e can lie anywhere in the range zero to 
one, zero corresponding to pure radial oscillations with vanishing net charge. Naively one 
might expect the charge of the resultant R-balls to be given by eR c . This is mainly true. 
However, once e drops below about 0.06, both positive and negatively charged R-balls are 
formed in comparable quantities and the typical charge becomes constant 
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4.2 Decay 



Bosonic decay modes do not exist when the vacuum of the low energy theory preserves 
R-symmetry (models which spontaneously break R-symmetry will be discussed shortly). 
The R-ball is already the lowest energy scalar field configuration for a given charge and 
gauge bosons have R-charge zero. However, several other decay modes do exist in most 
models; basically any light fermion with non-zero R-charge. Even so decay is slower than 
one might think due to the Pauli exclusion principle. Within an R-ball a Fermi pressure 
opposes the creation of fermions so, effectively, R-balls only evaporate from their surface. 

Consider first the couplings in the superpotential (12. 8p . In the absence of fine tuning one 
expects all particles to have tree level masses of order /x, or zero if imposed by symmetry. 
Since the condensate field mass m is generated at one loop it is comparatively suppressed 
and only decays to the massless fermions are kinematically allowed. Any massless fermion 
must be a null eigenvector of the fermionic mass matrix, given by \i in an R-symmetric 
model when (X) = 0. Consequently its scalar partner is massless (see ref. [I] or section 
[2j and the superfield must be a null eigenvector of A as well. Therefore massless fermions 
are forbidden from coupling to the condensate field at tree level. The condensate can only 
decay to the fermionic components of the <^'s at tree level if R-symmetry is spontaneously 
broken. 

A second option for R-ball decay is the gravitino. The goldstino has R-charge +1 
and is originally exactly massless, but gets eaten by the gravitino to pick up a mass 
m 3/2 ~ fr/Mp. Decay then proceeds via goldstino interaction terms [5THI5T] 

2 

£ D -^X^GG (4.13) 

JT 

with G denoting the gravitino. As the condensate field mass is generated at one loop 
the coupling is effectively loop suppressed. Loop decays to other light fermions (that 
may be external to the hidden sector) should therefore be considered too. One may also 
expect R-violating modes induced either by gravity, which is not excepted to respect any 
global symmetries, or the microscopic theory, which may break R-symmetry. These decay 
channels are suppressed by at least the cutoff scale so are not expected to be important. 

Rigorously estimating the lifetime of an R-ball in the thick wall limit is a difficult 
task [52] so here we shall take a simpler approach. The thin wall limit [SB] is not appropriate 
for the R-ball solution given in (14 .6p but the principles applied to bosonic decay modes in 
ref. [54] are. First consider tree level decays such as those to gravitinos. For a fermion 
coupling directly to the condensate field with strength g (equal to m 2 j fx for gravitinos) 
the penetration width is 1/gX. The fermion is produced inside the R-ball at some radius 
r so one must have 

— >r R -r (4.14) 

if it is to escape. Otherwise the Fermi pressure prevents decay. Since we are dealing 
with the interior of the R-ball, the large X component of eq. (14. 6 p is appropriate and the 
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inequality becomes 

~ u R r u R 1 

sm cuijr < — — ~ — = (4-15) 

ojrt gX R gv /R 

using urTr ~ 7r and, for the last term, eq. (14. SJ) - When the charge is large (specifically 
i? 3> g~ 2 ) this inequality can only satisfied around the surface of the R-ball, r = r\R. One 
thus finds 

R^-4ntu R drr 2 TX 2 (4.16) 

Jr R 

where T is the decay rate at a given point inside the R-ball and oorX 2 is the local charge 
density. 

Decreasing the charge of the R-ball by two (a single condensate quanta) liberates energy 
of order ur, yet fermions coupling to the condensate field at tree level gain local masses of 
order gX. Unless gX < ojr a local, tree level decay is kinematically forbidden and decay 
can only occur via heavy particle loops. On the surface X = a/7t by definition so the 
inequality is already satisfied there as long as g 2 < \R~ 1 ^ 2 . If not, the crossover occurs at 
a radius 



r-B for 7 ~l + _l„^_j (4.17) 

where the expressions for the small X component of eq. (14.61) and the R-ball parameters 
(14. 8 p have been utilised. Large R-balls have tr m tt/ujr ^> 1/m so the second term 
is subdominant and we can approximate 7 by one, meaning that gX < ujr everywhere 
from the surface outwards. The decay rate is therefore constant at its tree level value 
T « g 2 ujR/ 4tt. Essentially the R-ball solution decays very rapidly beyond the surface, so if 
gX < ujr is not already satisfied there it soon will be. Putting it all together we find 

R ~ (l + 2mTR + 2m 2 r 2 R ) ~ (4.18) 

4m 3 K RJ 2m K J 

where the last equality again follows from the fact that mrR m ixmjbjR ^> 1 for large 
charge^). 

Loop decays tell a similar story. As long as the mass of the heavy particle propagating 
around the loop is greater than the mass it acquires from its tree level coupling to the con- 
densate, the result T m g 2 uiR/ATi is unchanged. One simply replaces g with the appropriate 
effective coupling. Fortunately this will always be the case just beyond the surface for the 
reasons mentioned above. Coupling fields in the loop to light fermions with strength h 
thus gives a decay rate T ps g 2 h 4 u R /16iT 2 fj J 2 and, subsequently 

R y - — IILr^ 1 / 2 , 4.19 

8m/i 2 v ' 



8 In more detail, the decay rate starts off constant but, as the charge of the R-ball decreases, the other 
two terms grow and decay speeds up. Consequently R-balls decay more quickly with decreasing size, as 
was seen numerically for Q-balls in ref. |52j . 
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When T 2 < / the tree level decay rate is easily integrated to find an approximate 
lifetime 

2m A 4 

using eq. (I4.10p to set the initial charge and eq. (I2.10p to eliminate m. If loop decays 
dominate this result becomes 

16m/i 2 3/2 A 6 /i 
RS ~ 37rA 3 /i 4 /l ~ 10 3 AV^ 4 / 4 ' 

Otherwise the lifetime depends on the thermal history of the universe. Decay is quicker and 
the initial charge is smaller when the temperature is greater than yff, due to eq. (I3.16p . 
A lower bound is thus found by fixing fr at Tj, whereupon one can simply replace / with 
T 2 in the above expressions. 

4.3 Low energy theories without R-symmetry 

In O'Raifeartaigh models that spontaneously break their R-symmetry the initial situation 
remains broadly similar. Immediately after fragmentation the condensate moves in the 
logarithmic regime of the effective potential (I2.13p . This conserves R-charge so supports 
the formation of R-balls as before. The main difference is that there may now be R- violating 
bosonic decay modes available, which are not stifled by a Fermi pressure so can shorten 
the lifetime. Expressions for the fermionic decay rate are unaffected (we assumed nothing 
about charge conservation in the previous section short of establishing what decay modes 
were allowed) and the effect of bosonic modes can be estimated in a similar manner [53] . 

The key difference is that the lower limit of the integral (I4.16P is zero in the absence of 
Fermi pressure, so tree level decays acquire an extra contribution from the interior of the 
R-ball 

rr R / G 2 2 pr R rr R \ 

ARpa-Aiiu R J drr 2 TX 2 w - g W R \^—^J drr 2 + J drr 2 X 2 J . (4.22) 

Between r = and ^tr decay is at one loop via particles of mass gX, giving T m 
g 4 u R /167T 2 X 2 , whereas between and the surface we instead have V pa g 2 uj R /Ait. 
Due to the arguments associated with eq. (14.171) one expects 7 < 1, with equality when 
g 2 > \R~ l l 2 . For g 2 < XR" 1 ' 2 the crossover occurs inside the R-ball where the solution 
changes more slowly so one expects the deviation from 7 w 1 to be significant. Substituting 
in the large X component of eq. (14.61) and r R pa tt/ur the increase is given by 

AR ~ 1 y - - - — ^ 2vr 1 - 7 + sin 2tt7 

12 Aujr v v ,; /; 

~ - 9 Y^ r /? 3/4 (2 7 r(l - 7) + sin27T7). (4.23) 
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Large R means small ojr (specifically ur <C m, implying that cjr <C fr/i 71 ) whereupon 
the first term on the first line can be neglected. The second term, however, is extremely 
significant for large R-balls as it dominates the overall decay rate. 
Loop decays to bosons are similarly enhanced, contributing 

AR ~ - 92h4 i [ XfT f /2 R 1/4 (2ti(1 - 7) + sin2vr 7 ). (4.24) 

The positive exponent attached to R means that even one loop bosonic channels can be 
more important than tree level fermionic ones if the charge is large enough. However, note 
that if g 2 > XR" 1 / 2 (or g > fi/y/Jr for loop decays) 7 ~ 1 as before. Decay to bosons is 
then suppressed everywhere inside the R-ball due to the mass they acquire from condensate 
field couplings. Evaporation is only from the surface as per the R-symmetric case and there 
is no longer a significant enhancement. 

Assuming that g 2 < \R~ l l 2 and T 2 < f the revised lifetime for tree level decay modes 
is easily evaluated to be 

trb ~ ~A^R 1/A ~ T^7 ( 4 - 25 ) 



using eq. (j4.10p to set the initial charge and assuming that (27r(l — 7) + sin 2^7) /4 ~ 1. 
For loop decays one finds 

TRB 64 ^ . R^ „ (426) 

RB 3A 2 /i 4 (A/) 3 / 2 \ 3 g 2 h*P 1 J 

Just as in the R-symmetric case one can find a lower bound on the lifetime when the 
formation temperature is greater than y/J by replacing / with T 2 in these expressions. 
Of course, all of this assumes the existence of light bosons that the condensate is able to 
decay. It could be that all bosons are too heavy. Or the model parameters could conspire 
to enable the fermionic decay rate to overtake the bosonic one. In either case the lifetime 
reverts to that given in the previous section. 



5 R-ball phenomenology 

The cosmological behaviour of R-balls is wide ranging. To give a general overview of 
the features various models can exhibit, this section will focus mainly on O'Raifeartaigh 
models obeying the following criteria. The strength of all couplings in the superpotential 
( 12. 8 p will be taken to be of order A, whereas 'naturalness' suggests that the two scales in 
the superpotential ( 12. 8 p should be similar, i.e. /i 2 ~ fr, accounting for thermal corrections 
with eq. (I3.16p ^i We will further assume that all degrees of freedom in the hidden sector 



are heavy with masses /x, other than the condensate field. 

R-ball formation takes place in O'Raifeartaigh models whenever the condensate field 
(the scalar partner of the Goldstino) picks up a tachyonic soft mass due to couplings 



'increasing \i above this scale increases the tree level decay rate but decreases the one loop decay rate. 
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to the inflaton. When it does the characteristic R-ball scale A, defined in eq. ( 13. 7J) , is 
parametrically between the cutoff scale of the model and the SUSY breaking scale. The 
temperature of the universe at formation (I3.18P and typical charge ( I4.10p of large R-balls 
are then given in the following table 



Domain 


D l 


D 2 


D 3 


Tf 


/A/AfpN 1 / 2 
\ 4ttA / 


/A/A/pT^ 1 / 4 
\^ 4?rA J 


\ 4?rA J 


R 


A 4 


A 4 


A 6 


10 5 A 2 / 2 


10 5 A 2 / 2 


10 3 A 4 A/ p T 4 



(5.1; 



Domain boundaries are determined by when R-ball formation takes place relative to the 
decay of the inflation, and whether thermal effects dominate the effective potential during 
this process. They are defined by 



Di : T f <T h ,Tf<f 



D 2 : T f >T h ,T f 2 <f 



Ih ■■ // > V/. . Tj ■ J 



A > max 



XfMp XM P ffMp 



^1 ' 



XfM P 



> A > max 



1/3 

1 V 4ttA ') 
\MpTl (fMpTl 



mm 



\MpTl XM P 



Airf 



Air 



4nf ' 
> A > 



4vrA 3 
4vrA 



1/5' 



1/3 



(5.2) 



where Th denotes the reheat temperature of the universe after inflaton decay. Outside of 
these domains the condensate is eroded before R-balls are able to form. For gravity medi- 
ated SUSY breaking the reheat temperature must lie in the range 10 MeV < Th < 10 6 m 3 / 2 
whereas for gauge mediation the bound is 10 MeV < T h < 10 14 (m 3 / 2 /GeV) _2 GeV. 
The initial energy density stored in R-balls ( 14. lip is given by 



10A 2 



(5.3) 



(or one if A ~ Mp) and scales with the subsequent expansion of the universe (I4.12p as 

{10A 2 /MJ, for Tf > T > T h 

(T h /T)(10A 2 /M P ) for T h > T > T e (5.4) 
(T,/T e )(10A 2 /M 2 ) for T e > T 

where T e ~ 1 eV is the usual temperature at which the universe becomes matter dominated. 
If Tf < T h one replaces T h with Tf in the second equation and, for A G D 3 , the density 
is multiplied by ^J~fr|Tf due to the temperature dependence (13.161) of the R-ball energy 
(USD- 

From here on the analysis becomes strongly dependent on the particulars of the model 
and the available decay modes. These in turn are specific to the mechanism chosen to me- 
diate SUSY breaking. The main examples we will consider are gravity and gauge mediated 
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SUSY breaking, for preserved and spontaneously broken R-symmetry, both for strongly 
and weakly coupled superpotentials. Some general features that crop in these models are 
as follows. 

Perhaps of most immediate interest is the idea that R-balls are long lived and still exist 
today, contributing to the dark matter density of the universe. According to eqs. (I4.20p . 
(14.2 ip . ( I4.25P and (I4.26P long lifetimes correspond to a high cutoff scale (or an exact R- 
symmetry), a small SUSY breaking scale (or equivalently a small gravitino mass) and/or 
a weakly coupled condensate field. If sufficiently long lived, non-relativistic R-balls behave 
as cold dark matter with a density given by eq. (15. 4p . One must find f2^(T e ) < 0.22 if 
they are not to exceed the total observed value. Note that light gravitinos, one of the 
requirements for R-ball dark matter, are specific to gauge mediated SUSY breaking, hence 
one does not expect this kind of dark matter in models of gravity mediation. 

If not sufficiently long lived to survive until the present day R-balls must, obviously, 
decay. The effect this has on the visible sector varies and could in principle result in either 
heating or cooling of the universe. Each charge two quanta of condensate that decays 
carries away energy 



using eqs. (14.81) and (I4.10p . Regardless of what the decay products actually are, this sets 
their maximum mass and characteristic temperature. 

Since R-balls are not in thermal equilibrium with the rest of the universe their decay 
products will initially be out of equilibrium as well. Large R-ball density (15.41) during decay 
thus results in a secondary reheating, or cooling of the universe. An immediate consequence 
is that, if one occurs, an R-ball dominated epoch decouples the generation of visible sector 
matter and radiation from the dynamics of the inflaton. R-ball, rather than inflaton, decay 
can be responsible for the present contents of the universe, with the new reheat temperature 
AE obeying the constraints applied to the original one TV If, on the other hand, the R- 
ball density is small the temperature of the universe is unchanged. The decay products 
are brought into equilibrium with everything else, unless they are incapable of maintaining 
thermal equilibrium in which case they simply update the relevant relic abundance by a 
small amount. Either way the scenario is of limited interest. 

R-symmetric models permit decays to light fermions alone. The only such fermion 
that can couple to the condensate field at tree level is the gravitino but, owing to the 
loop suppressed mass of the condensate field, decays at one loop may also be important. 
This is particularly so in the absence of light hidden sector states whereupon these decays 
are to visible sector gauginos (which have R-charge +1). When the O'Raifeartaigh model 
spontaneously breaks R-symmetry R-balls can also decay at one loop to visible sector 
gauge bosons, sfermions or even fermions. Alternatively, if light hidden sector degrees of 
freedom are included, one generally finds an increased decay rate as R-balls can decay 
directly to said fields rather than via loops or gravitinos. However, the observable effects 
in this scenario are highly model specific and depend on the exact properties of these new 




for A e D U D 2 
for A e D 3 



(5.5) 
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degrees of freedom so will not be discussed here. 



5.1 Gravity mediation 

In pure gravity mediation the O'Raifeartaigh model communicates with the visible sector 
only through gravitational interactions. Gravitinos couple to the condensate at tree level 
with strength g ~ m 2 / _/r ~ A 4 /167r 2 so lead to an R-ball lifetime (14.20)) 



10A 4 /A 8 (m 3/2 M P ) 5 / 2 for A e D x , D 2 
10A7A 8 T| for AeD 3 ^ 



in R-symmetric vacua, where this is the only available decay mode. The only light states 
accessible to the condensate at one loop are those in the visible sector and it must be 
gravitinos propagating around the loop. Loop decay rates are therefore negligible as they 
contain higher orders of the already small gravitino couplings. 

Spontaneously broken R-symmetry opens up the possibility of decay to visible sector 
gauge bosons via gravitino loops. This decay rate is only significant if the R-ball interior 
contributes, i.e. g < 1 (see section |4T3|) . which is always true for the small gravitino coupling. 
Gravitino couplings to gauge bosons go like h ~ m 3 / 2 /Mp leading to the shortened R-ball 
lifetime ( jjjZg]) 

f 10 5 A 3 M 2 /A n m6 forAGDi,/^ 
trr ~ < 1 (5 7) 

1 10 5 A 3 Tf/X u ml /2 iorAEDs V ; 

provided this expression evaluates to less than (I5.8p . otherwise decay remains to gravitinos. 

We can use these lifetimes to slice up the O'Raifeartaigh model parameter space (m 3 / 2 , 
A, A) into different regions of interest. It actually turns out that R-balls in gravity mediated 
SUSY breaking do not have a wide range of phenomenological consequences; they can only 
result in secondary reheating of the universe by gravitinos. This is because loop couplings 
remain sufficiently suppressed over the entire parameter space to favour gravitino decays, 
and the high SUSY breaking scale tends to keep their lifetime short. Weak coupling does 
lead to a small region supporting long lived R-balls, but the density there is too high to 
be compatible with the observed cold dark matter density. Some R-symmetric examples 
are given in figure I5TT1 R-breaking O'Raifeartaigh models still result in R-ball formation, 
but their density is typically too low to have any interesting consequences. 



5.2 Gauge mediation 

A combination of messenger loops and a low gravitino mass in gauge mediation lead to a 
much more interesting phenomenology. Although gravitinos are still the only light fermion 
coupling to the condensate field at tree level, there are now one loop couplings via messen- 
gers to visible sector gauge fields. These loops are much more significant than gravitino 
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RS-GR, T h = 10 6 m 3/2 , A = 1 RS-GR, T h = 10 MeV, X - 10" 3 




2 3 4 2 3 4 



Figure 5.1: R-balls in R-symmetric, gravity mediated SUSY breaking as a function of 
log (m3/2/GeV) (horizontal axis) and log(A/GeV) (vertical axis) for various values of A 
and Th- In region 1 they have no significant effects. In region 2 they reheat the universe 
by decaying to gravitinos. Outside the shaded region R-balls do not form due to thermal 
effects, or they contribute more than the observed cold dark matter density. 



loops as messengers couple to the condensate field with strength g ~ A, then through vis- 
ible sector gauge couplings with strength h ~ 1. If gauginos are kinematically accessible, 
i.e. AE > 1 TeV ( 15.51) . and decay to them is faster than to gravitinos, the R-ball lifetime 
( ECHD is 

f A 6 /10 3 A 4 (m 3/2 M P ) 7 / 2 for AeD u D 2 
r R s ~ < (5.8) 
1 A 6 /10 3 A 4 T; for A 6 Z) 3 

for an R-symmetric O'Raifeartaigh model. Otherwise the lifetime is determined by the 
decay rate to gravitinos as in eq. (15. 6p . 

R-breaking O'Raifeartaigh models are even more sensitive to loop decays. Now visible 
sector gauge bosons are always accessible and dominate decay if either eq. (15. 8p or 

f 10A 3 /A 5 m2 M 2 forAe£>i,£> 2 
trb ~ < ' (5.9) 

) 10A 3 /A 5 T; for AG D 3 . 1 ' 

(using eq. ( I4.26P ) is less than eq. (15. 6p . In either case R-balls preferentially evaporate to 
gauge bosons rather than fermions due to the extra contribution from the interior. 

The parameter space (m 3 / 2 , A, A) can be split as for gravity mediation but we now 
find more variation in the results. R-symmetric O'Raifeartaigh models demonstrate both 
gravitino reheating and cooling, gaugino reheating and R-ball dark matter. Dark matter 
corresponds to a small original reheat temperature so as to minimise the boost given to the 
density (15. 4p . Otherwise the regions of parameter space coincident with long lived R-balls 
tend to be associated with too high a density. Furthermore the gravitino mass is small in 
these regions so as to maximise the lifetime. 
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RS-GA, T h = 10 14 (m 3/ 2/GeV)" 2 GeV, A = 1 RS -GA, T h = 10 14 (m3 /2 /GeV)" 2 GeV, A = 10 " 3 RS-GA, T H = 10 MeV, A = 10" 3 
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Figure 5.2: R-balls in R-symmetric, gauge mediated SUSY breaking as a function of 
log (m3/2/GeV) (horizontal axis) and log(A/GeV) (vertical axis) for various values of A 
and Th- In region 1 they have no significant effects. In region 2/3 they reheat/cool the 
universe by decaying to gravitinos. Region 4 corresponds to reheating by gauginos and 
region 5 to R-ball dark matter. Outside the shaded region R-balls do not form due to 
thermal effects, contribute more than the observed cold dark matter density or their decay 
violates the bounds on the reheat temperature. 



For most of the range in which R-balls live long enough the gravitino mass is less than 
a few keV. Its contribution to the overall dark matter density is therefore small [S3], but 
R-balls can easily account for the entire Qdm ~ 0.22. If they do (the top of region 5 in 
figure 15.21) their charge, size and energy (14.81) are in the ranges 

10 29 < R < 10 49 , lO" 15 m < r R < 1(T 4 m, 10 28 GeV < E R < 10 38 GeV. (5.10) 

Individual R-balls can therefore be anything from fermi to micro scale and, regardless of 
size, are very dense objects. Outside of the dark matter region gravitino decay is the 
norm for R-symmetric models due to the kinematic constraints on forming gauginos, but 
there are regions of parameter space favouring evaporation to gauginos. Three examples 
of R-balls in R-symmetric, gauge mediated SUSY breaking are provided in figure 15.21 

O'Raifeartaigh models with spontaneously broken R-symmetry have a somewhat dif- 
ferent phenomenology. In all regions where R-ball decay has a significant impact on the 
evolution of the universe decay is to visible sector gauge bosons, and can result in ei- 
ther heating or cooling. Cooling takes place if the original reheat temperature is high 
and A < 10 18 GeV (the decreased lifetime wins out against the increased energy), whereas 
heating occurs elsewhere. There also remains a small region supporting R-ball dark matter, 
with parameters skewed towards the large R-balls of eq. (I5.10p . Some example of R-balls 
in R-breaking, gauge mediated SUSY breaking are shown in figure 15.31 



23 



RB-GA, T h = 10 14 (m 3/2 /GeV)- 2 GeV, X = 1 RB-GA, T h = 10 MeV, A = 10" 3 




-4 -2 2 4 -12 -8 -4 4 



Figure 5.3: R-balls in R-breaking, gauge mediated SUSY breaking as a function of 
log (m3/2/GeV) (horizontal axis) and log(A/GeV) (vertical axis) for various values of A 
and Th- In region 1 they have no significant effects. Region 5 corresponds to R-ball dark 
matter. In region 6/7 they reheat/cool the universe by decaying to visible sector gauge 
bosons. Outside the shaded region R-balls do not form due to thermal effects, contribute 
more than the observed cold dark matter density or their decay violates the bounds on the 
reheat temperature. 



5.3 Detecting R-balls 

Experimentally, dark matter R-balls would be challenging to observe. They cannot be 
produced in colliders, but one might hope for a dominant decay mode to visible sector 
particles that can be observed in some other type of experiment. Gravitino decays and any 
decays to the hidden sector are thus ruled out, but decays at one loop via the messengers of 
gauge mediation are a possible candidate. If R-symmetry is preserved the energy released 
in the decay of a single condensate quanta (15.51) must be sufficient to produce a pair of 
gauginos, but if R-symmetry is spontaneously broken R-balls can always decay to massless 
visible sector gauge bosons. From eq. (I5.10p the expected energy range for R-ball dark 
matter is 

10" 11 GeV < AE < KT 1 GeV (5.11) 

which is well below the gaugino mass. Ergo dark matter R-balls can only be observed 
through decay to visible sector particles if R-symmetry is spontaneously broken. 

Q-ball detection techniques, as studied in refs. |15 ]I5^Hd^] , do not apply to R-balls as the 
condensate field is not charged under the visible sector gauge group. One could perhaps 
search for them using direct detection experiments though. Visible sector matter will 
scatter elastically off R-balls at one loop, through penguin diagrams containing messengers, 
for example. Owing to the classical nature of R-balls this process is likely to be somewhat 
non-standard and could produce a distinctive signature. The details are left for future 
work. 

Evidence for R-balls decaying before the present day would surely necessitate the incep- 
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tion of a test with a more cosmological nature. Their formation and decay are potentially 
significant events in the evolution of the universe so may well have left an imprint on some 
large scale, cosmological observable. However, the details of such a test are beyond the 
scope of this work. 

6 Summary 

Condensates forming along flat directions of O'Raifeartaigh models can have a significant 
impact on the evolution of the universe. They are somewhat generic, emerging in any model 
where the flat direction acquires a tachyonic soft mass through couplings to the inflaton. 
When a condensate does form it eventually fragments into non-topological solitons with 
conserved R-charge, known as R-balls. These objects are large, classical configurations and 
allow an approximate, analytical description. Formation is insensitive to whether or not 
R-symmetry is spontaneously broken, but decay is not. 

Depending on the scale up to which the O'Raifeartaigh model is valid, the scale of 
SUSY breaking and the strength of tree level couplings in the superpotential R-balls result 
in a variety of phenomena (figures 15.21 and 15. 3p . In gravity mediated SUSY breaking 
they can reheat the universe through gravitino decays. In gauge mediation they provide 
a good dark matter candidate if stable, or decay to gravitinos, gauginos or visible sector 
gauge bosons, either reheating or cooling the universe. Both mediation mechanisms enable 
one to decouple the generation of visible sector matter from inflaton dynamics, instead 
using R-balls to reheat the universe. 
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